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Abstract 



A pentagon in the plane with fixed side-lengths has a two-dimensional shape space. 
Considering the pentagon as a mechanical system with point masses at the corners we 
answer the question of how much the pentagon can rotate with zero angular momen- 
tum. We show that the shape space of the equilateral pentagon has genus 4 and find 
a fundamental region by discrete symmetry reduction with respect to symmetry group 
Dq. The amount of rotation A9 for a loop in shape space at zero angular momentum 
is interpreted as a geometric phase and is obtained as an integral of a function B over 
the region of shape space enclosed by the loop. With a simple variational argument we 
^ determine locally optimal loops as the zero contours of the function B. The resulting 

shape change is represented as a Fourier series, and the global maximum of AO ~ 45° is 
found for a loop around the regular pentagram. We also show that restricting allowed 
^ shapes to convex pentagons the optimal loop is the boundary of the convex region and 

^ gives A6I f« 19°. 

(N 

• 1 Introduction 

^-H The possibility of achieving overall rotation at zero total angular momentum in an iso- 

lated mechanical system is surprising. It is possible for non-rigid bodies, in particular for 
^ systems of coupled rigid bodies, to change their orientations without an external torque 

^ using only internal forces, thus preserving the total angular momentum. The classical 

example of this phenomenon is exhibited by falling cats: a cat dropped upside-down 
without angular momentum will reorient itself by changing its shape and land on its 
feet, with roughly the same final and initial shape. The first theoretical explanation was 
given by Kane and Scher [6] , also see [lUl E] . 

Here we present a study of the equilateral pentagon in the plane where we permit 
the angles ipi, i = 1, . . . ,5 of adjacent edges to continuously change while the lengths of 
the edges are all fixed to 1. The reason the pentagon is chosen among other polygons is 
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because it has a two-dimensional shape space. By contrast, an equilateral triangle has a 
fixed shape; an equilateral quadrilateral can change from a square, through a rhombus, 
to a degenerate shape of a line. Adding one more degree of freedom makes it possible for 
the equilateral pentagon to achieve overall rotation at zero angular momentum through a 
periodic shape change. In general the side- lengths of the pentagon could be considered as 
parameters, but we restrict ourselves to the equilateral case, which gives some additional 
simplification and beauty through its discrete symmetry. 

Changes in size are irrelevant to our problem, so we use the word shape in the sense 
of congruence: two pentagons have the same shape if one can be transformed into the 
other by isometrics of the Euclidean plane, that is a combination of rotation, translation 
and reflections, generating E{2). Sometimes we will consider direct isometrics SE{2) 
only, omitting the reflections. Thus the equilateral pentagon is a mechanical system 
with symmetry. In this setting the overall rotation at angular momentum zero appears 
as a geometric phase. Symmetry reduction splits the dynamics into a motion in the 
symmetry group (translations and rotations) and a reduced system. The motion in 
the reduced system drives the motion in the group direction, and it is possible (in fact 
typical) that traversing a closed loop in the reduced system does not lead to a closed 
orbit in the group, see for example ISil^. The motion along the group can be split into 
a geometric and a dynamic phase, where the geometric phase does not depend on the 
speed at which the loop in the reduced system is traversed. In our case the translation 
is removed by going to the centre of mass frame, and since we are in the plane only a 
single angle 6 is needed to describe the orientation. Although the inspiration for this 
work was taken from the general modern theory of geometric phase [Si [71 [1], here we 
take an approach that can be understood with a minimal background in mechanics. 

The equilateral pentagon as a symmetry reduced mechanical system would perform a 
certain motion in (the cotangent bundle of) shape space if the angles were free to move. 
In many ways the system would then be similar to the 3- linkage studied in [2] , the four- 
bar linkage studied in [T2], or, e.g., the planar skater studied in [S]. Here, however, we 
take the point of view that the angles can be completely controlled by us, for instance 
with a motor per joint. The only constraint imposed is that the motion must be such 
that the total angular momentum remains constant (at value 0, in particular). A way 
to picture this is to think of the pentagon as a space station that has five motors at the 
joints and is floating in space without angular momentum. Controlling the motors, we 
are free to prescribe any motion in shape space. We can then ask "What is the optimal 
periodic shape change of the equilateral pentagon so that the overall rotation after 
one traversal of the loop in shape space is as large as possible?" Specifically we seek the 
global maximum of on the space of all finite smooth contractible loops. 

The plan of the paper is as follows: 

• Section 2 - We describe the equilateral pentagon and its shape space; 

• Section 3 - Using reduction by the discrete symmetry group -D5, we obtain a 
fundamental region of shape space and show how all of shape space is tiled by this 
fundamental piece; 
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• Section 4 - Explicit formulas for the moment of inertia, angular momentum and 
the rate of change of orientation are derived; 

• Section 5 - The geometric phase A9 is defined in terms of a line integral which is 
then converted into an area integral over the enclosed region of a scalar function 
B on shape space; 

• Section 6 - Wc show that the zero level of the function B gives the optimal loop, 
and we obtain a representation of the corresponding shape change in terms of a 
Fourier series; 

• Section 7 - Restricting to convex pentagons we show that the optimal loop for this 
sub-family is given by the boundary of the region of convex pentagons; 

2 Equilateral Pentagons 

The equilateral pentagon in this study has the following attributes: 

• Vertices are treated as point particles, each of unit mass; 

• Each edge is masslcss and fixed at unit length; 

• The angles between adjacent edges are allowed to change freely; 

Note that the family of all equilateral pentagons includes degenerate pentagons (e.g. 
an equilateral unit triangle with extra folded edges or a trapezium with one of the 
pentagonal angles taking vr ) and non-simple pentagons (e.g. the pentagram). 

We denote the vertices of the pentagon by G C where z G Z5, that is, the vertex 
indices are always modulo 5. We will represent the elements of Z5 by {1, 2, 3, 4, 5}, start- 
ing with 1. The standard colour code that we use for the vertices are {zi, Z2, 2:3, 2:4, 25} = 
{green, red, blue, black, yellow}, this colour code will also be adopted for the relative 
angles ■^j. The oriented edges of the pentagon are the vectors z^+i — Zi. 

A polygon is called simple if the edges do not intersect except at the vertices. The 
internal angle sum of a simple, n-sided polygon is (n — 2)7r, since its interior can be 
tiled by n — 2 disjoint triangles. This construction does not work when a polygon is 
self-intersecting. In fact, even the notion of an internal angle is not well defined in a self- 
intersecting polygon. For this reason we adopt a convention for measuring the angles of 
the equilateral pentagon that gives the natural result of it^ and ib|, correspondingly, 
for all of the relative angles of the regular convex pentagons and regular pentagrams, 
where the sign denotes orientation. We define the relative angle ipi as the amount of 
rotation needed to turn the oriented edge Zj+i — Zi into the negative of the previous 
oriented edge zi — Zi-i, hence e*'^»(zj+i — Zi) = —{zi — Zi-i) where i € Z5. As a result of 
this definition the relation between successive vertices G C is 



Zi+i = Zi-e '"^'{zi - Zj-i). 



(1) 
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For definiteness we use the principal argument Arg to define ipi from the vertices Zi, such 
that 



Arg 



Zi+l 



Z'i-l 



(2) 



The relative angles ipi define the shape of the pentagon. The shape of the pentagon 
is invariant under rotations Zi — )• e*^Zj. The absolute angle 9 is introduced to measure 
the orientation of the pentagon. We define 6 to be the angle between the positive x-axis 
and the edge Z2 — z\ measured counter-clockwise from the x-axis, such that Z2 = z\ + e*^. 

The shape of the pentagon is invariant under translations Zi — >• Zi^z. Initially we let 
the first vertex be arbitrarily located at z, but this z will later be eliminated by fixing 
the centre of mass at the origin. The vertices are: 



z\ 

Z2 

z-i 

Z^ 



z 

Z2 - e'^^e-'"^^ 

Z3 + e*^e-*(^2+V3) 



24 



g«Sg-«('/'2+l/'3+V'4) 



By a translation we can achieve zi + Z2 + z^ + z^ + z^ = Q so that the centre of mass is 
at the origin, and thus eliminating z gives: 



1 



zi 



Z2 



zz 



Zi 



Z5 



JO 



-4 + 3e" 



-iip2 



2g-*(^2+V'3) _|_ g-«(V'2+V'3+V'4)^ 



I _|_ 3g-*'/'2 _ 2e~*('^2+'/'3) _|_ g-«('/'2+'/'3+V'4)^ 



2g-#2 _|_ ^g-i{^2+^3) _|_ g-«('/'2+V'3+'04) 



I _ 2g-*'/'2 _|_ 3g-«(V'2+V'3) _ 4g-«(V'2+l/'3+V'4)^ 



(3a) 



(3b) 



(3c) 



(3d) 



(3e) 



In the following the symbol Zi will refer to these formulas. 

The ordered set of vertices Zi defined by ([s]) gives the equilateral pentagon modulo 
translations. The ordered set of relative angles ipi gives the equilateral pentagon modulo 
the special Euclidean group SE{2) of orientation preserving rotations and translations. 
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Thus the relative angles describe the labelled (and hence oriented) shape of the equi- 
lateral pentagon obtained by reducing the continuous symmetry SE{2). Later we will 
consider additional discrete symmetries to reduce further. They are the symmetry group 
of the equilateral (and equal masses) pentagon, and the reflection Z2 = E{2) / SE{2). 
The quotient by the full symmetry group x E{2) gives (unlabelled) shape up to 
congruence. 



2.1 Constraints 

When the vertices are placed ak, Zi = i and consecutive edges are connected with joints 
at Z2, 23 and 24 we obtain the so called 4-segment open linkage system. The shape space 
of this system is topologically a 3-dimensional torus, with angles V'S) ip4,. The 
open linkage can be closed by requiring that \zi — z^\ = 1, introducing the fifth edge of 
length 1 connecting 25 and zi. The closure constraint reduces the dimension of shape 
space to 2, and turns its topology into that of a rather complicated surface of genus 4 

laE]. 

The relative angles and Vi do not enter the equations for Zj, see ([3|, and are 
completely determined by the other angles. The equation for ip^ is obtained from ([T]) 
for z = 5, both sides are multiplied by e*('^2+V'3+'/'4)^ then the complex logarithm is 
taken to obtain 

V's = - Arg {l - e'^^ + e*('^3+V'4) _ gi(^2+V3+^4)^ (4) 

A similar calculation is used to show that 

V'l = Arg {l - + g-iW2+V3) _ e-»{^2+V3+V'4)^ _ (5) 

The above angles are all related by 

V'l + V'2 + V-s + ^^4 + V'5 = (1 + 2A;)7r, A: G Z . (6) 

The sum of all five relative angles ^ V'i takes values — Svr, — vr, vr, Svr, corresponding to 
k = —2,-1,0,1. The extremal values are achieved in the region around the regular 
simple pentagon where the pentagon remains convex, with either positive or negative 
orientation. Beyond these regions the pentagon may or may not be simple, but its angle 
sum remains constant as long as no additional stretched edge-configuration with ■0j = ibvr 
appears. We could define the orientation of a pentagon by the sign of ^ Vi- For shapes 
with ij^i = zbvr this orientation is undefined as nearby shapes have either sign. 

The complicated topology of shape space of equilateral pentagons arises because the 
relative angles ip2i V'3) ^'^d ip4^ are not independent. Given two arbitrary angles "02 and 
ip2,i the relative orientation of three of the edges of the pentagon is fixed as shown in 
Figure l|a) For the pentagon to be able to close with the two remaining segments, it 



is necessary that < \zi — zi\ < 2. There are four possible cases when attempting to 
determine ^^4 from -02 and -03: 

• 1^4 — zil = 0: ^4 is undetermined, there are infinitely many solutions. 



2 EQUILATERAL PENTAGONS 



6 



Z5 




(a) 3 segments of the pentagon shown. (b) With -02 and xps specified, there are two possible 

pentagons in general. 

Figure 1: Pentagon construction by specifying two successive relative angles. 



• < |z4 — zi| < 2: the generic case where there are exactly two solutions for -04 as 
illustrated in Figure []|[b)[ 

• 1^4 — zi| = 2: the special case that there is a unique solution for ^4. 

• 1^4 — zi| > 2: there is no solution since the pentagon cannot close. 

The existence of the first case with infinitely many solutions makes the shape space more 
complicated than some gluing of the 2-torus {ip2,ip3) of the 3-segment open linkage. 
Starting with the closure condition 

\Z5 - zi\^ - 1 = 0, (7) 

and inserting ^ gives the relation between the relative angles 1^2, ip3^ -04 as 

3 — 2 cos ^2 — 2 cos Tps — 2 cos ip4 + 2 cos {lp2 + tp3)+ 
+2 cos (V's + i^i) - 2 cos {lp2 + ip3 + ipi) = . 

This equation defines the shape space as a 2-dimensional sub-manifold of the 3-dimensional 
torus, the shape space of the 4-segment open linkage. Equation ([7]) can be rewrit- 
ten as XC + CC + X-^C = where X = e*'^*, C = -1 + e'^^ - e^^+'^-'\ such that 
CC = 1-24 — ^iP = 3 — 2 cos 'i/'2 — 2 COS V's + 2 COS {ip2 + V's)- Assuming |C| / , and using 
the polar form of ( gives 

ip4 = zbarccos ^""^^ ~ argC s.t. V4 £ {—'T^^M ■ (9) 

When \Q = \z4 — zi\ = the value of tp4 is undetermined. This occurs only when 
■02 = V'3 = ±7r/3. Thus, solving of ([T]) reflects exactly the four cases discussed above. 
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2.2 Shape Space 



Although the relative angles ipi are easy to visualise and interpret geometrically, 
algebraic equations can be simplified with the following affine transformation 



the 











. "3 . 





-1 








V'2 




" " 




V'3 


+ 


vr 




_ Ipi _ 




_ _ 



Ail^ + b. 



(10) 



The determinant of ^ is the negative sign means that the orientation is reversed, 
and the factor ^ means that the area is halved when we do the transformation. That 
means if we use the same natural domain (— vr, vr] for both the ^ and a-coordinates, then 
we have a double covering in the a-coordinates. That is, every possible pentagon in the 
■i/;-coordinates occurs exactly twice in the a-coordinates. Since the angles 'i/'i are defined 
modulo 2vr this induces an equivalence relation for the angles a^, which is given by 

(ai, a2, as) = (ai, a2, as) — ((z + 2j + /c)vr, (i-l- A;)vr, (i — A;)vr) for i,j,k£lj. (11) 

Using this equivalence relation the double covering can be removed by restricting the 
fundamental domain to our choice of as E [0,vr). 

The transformation is invertible, and its inverse is 



1p2 







-1 


-1 " 




ai 








-1 


1 







a2 — vr 


mod 27r 


_ tp4, _ 







-1 


1 




as 





(12) 



For convenience, both coordinate systems will be used. The geometrical interpreta- 
tions will be done in the ^-coordinate system while the algebraic calculations will be 
done in the a-coordinate system. In a-coordinates the equation for shape space ([s]) 
simplifies to the symmetric form 



C(ai, a2, as) := 3 + 4cosai cosa2 + 4cosai cos as -|- 4cosa2 cos as = 0. 



(13) 



The corresponding surface C(ai,a2,as) = is shown in Figure [2| It is reminiscent of 
the I-WP surface from minimal surface theory [5]. Our defining equation is like a low 
order Fourier approximation of the I-WP surface. 



From equation (13), as can be uniquely expressed in terms of ai and a2 as 



as 



arccos 



-3 — 4cosai cos a2 
4(cos ai + cos 02) 



(14) 



A unique solution is obtained for as as we are strictly in the domain [0,vr). Accordingly 
ai and a2 are good local coordinates almost everywhere on shape space. The exceptions 
are at the boundaries where as = 0, vr, and, in particular, when (ai, 02) takes the values 
(ibvr/6, ib5vr/6) and (±5vr/6, ibvr/6). Over these points there is a whole line in the surface 
corresponding to infinitely many solutions of as- This corresponds to the case when the 
fraction in (14) is of the form 0/0. All other points for which cosai +cosa2 = are not 
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Figure 2: Shape space in a-coordinate system. In ai and 02 there are periodic boundary 
conditions. The gluing from 03 = to 03 = vr is done after a shift by vr in ai and 02- 
This is a surface of genus 4. 

on shape space at all, so we do not need to worry about them. We will see that after 
discrete symmetry reduction it is not necessary to consider the singular lines anymore 
(where the exceptions for ai and 02 occur). 

Perhaps the simplest way to understand why the genus is 4 is to compute the Euler- 
characteristic from the picture taking the gluing into account. Alternatively we can 
compute it using Morse theory from the critical points of a smooth function defined on 
the surface, we will follow this approach with the moment of inertia later on. 

Notice that our definition of shape space is the set of all equilateral pentagons with 
labelled (or distinguishable) vertices. We use the colours green, red, blue, black, yellow 
to designate the vertices zi, Z2, z^, z^, Z5, respectively. In the standard notion of con- 
gruence of polygons the vertices (and sides) are unlabelled and hence indistinguishable. 
Considering the polygon as a mechanical system (with potentially different masses at 
the vertices and different moments of inertia for the edges) we obtain a description of the 
space of labelled equilateral pentagons. As it turns out, all quantities we are interested 
in here are given by functions on the unlabelled shape space, because the masses at 
the corners and the side-length are all equal. The reduction from labelled to unlabelled 
shapes is a discrete symmetry reduction, which we are going to study next. 

3 Discrete Symmetries 

Consider an arbitrary equilateral pentagon with 5 labelled vertices, e.g. distinguished by 
colours. The action of the group -D5 on the labelled shape space is generated by vertex 
rotations and vertex reflections. By a vertex rotation we mean a cyclic permutation 
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of the vertices, that is R{zi, Z2, Z3, Z4, z^) = {z2, Z3, Z4, Z5, zi). By a vertex reflection 
we mean a permutation that reverses the order of vertices and fixes a chosen vertex; 
choosing to fix Z3, we get V{zi, Z2, zs, Z/^, z^) = (z^, z^, z^, Z2, zi). Both operations leave 
the set of vertices Zi in the plane fixed; they merely permute the vertices in a way that 
preserves neighbours. Thus the unlabelled shape of the pentagon is fixed as well, but 
the orientation induced by the labelling is reversed under V. For a generic shape both 
operations do change the labelled shape, that is there is no element in E(2) that undoes 
the action. Special labelled shapes are fixed under subgroups of D^, for example the 
labelled regular pentagons are invariant under vertex rotations R, since up to a geometric 
rotation G SE[2) it is the same labelled shape as before vertex rotation. Symmetry 
reduction of the labelled shape space allows us to define a fundamental region in which 
each unlabelled shape is represented exactly once. Since the action of is not free, the 
unlabelled shape space is not a smooth manifold but just an orbifold, with singularities 
at the shapes that have higher symmetry, and hence non-trivial isotropy. 

If we write the relative angles as a vector = {ipi,ip2, ^3) ^'4) V's)*; then the symmetry 
operations can be represented by multiplication of ■0 by matrices R, for vertex rotation, 
and V, for vertex reflection (fixing vertex 3) where 
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The additional minus sign in V follows from the definition of ipi in ([2]). The group 
is generated by R and V with presentation {R,V\R^ = V'^ = id,VR = R^^V). The 
group -D5 acting on the relative angles ip leaves the unlabelled (and un-oriented) shape 
invariant. In general it changes the labelled shape; normally there are 10 (the order of 
D5) different labelled shapes corresponding to the same unlabelled un-oriented shape. 

There is another discrete symmetry because by using the relative angles ipi to describe 
the shape we have reduced by SE{2), but not by E(2). Hence reflections about a 
line through the origin give another discrete symmetry Z2 = E{2) / SE{2). We call it 
the mirror reflection symmetry M. The action on the vertices is M{zi, Z2, z^, Z/^, z^) = 
{zi, Z2, Z3, Z4, Z5), where the overbar denotes complex conjugation. As a matrix acting 
on the space of angles ip we simply have M = —id. Unlike R and V the operation M, 
in general, changes the set of vertices Zi, even modulo SE{2). 

Labelled polygons have an orientation induced by the labelling, while for unlabelled 
polygons an orientation may be kept track of by orienting the edges with an arrow. Con- 
sidering oriented pentagons, both M and V reverse the orientation. Simply forgetting 
the labels of the vertices gives an unlabelled un-oriented shape. This corresponds to 
reduction by the full group -D5 x E{2), which gives a shape in the classical sense that 
two shapes are the same if they are congruent (and have the same size). We will see that 
unlabelled oriented shapes up to congruence are obtained from reduction by a slightly 
different group D+ x SE{2). 
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Combining all three discrete symmetries gives x Z2 which is isomorphic to Dio 
with presentation {MR,V\{MRy^ = V'^ = id,V{MR) = {MR)-^V). A subgroup 
of Dio, different (but isomorphic) to generated by R and V , is obtained from the 
generators R and MV: D+ = {R,MV \ R^ = {MV)"^ = id, {MV)R = R-^{MV)). The 
superscript + indicates that the action of this group preserves orientation as both R 
and MV preserve orientation. The group D'^ is used for the symmetry reduction in 
the next chapter, and the resulting reduced shapes are unlabelled oriented shapes. In 
the following we will drop the qualification unlabelled and simply talk about oriented 
shapes. 



3.1 Discrete Symmetry Reduction 

We now construct a fundamental region of labelled shape space such that the whole 
surface is obtained as the orbit of this fundamental region. This fundamental region 
will contain every oriented shape exactly once. If necessary, a final reduction by V to 
remove the double covering of un-oriented shapes from this region may be performed to 
halve the fundamental region. 

The important objects in symmetry reduction are isotropy subgroups. The isotropy 
group of a point x in shape space is defined as Gx = {g ^ G : gx = x}. A shape has 
high symmetry if it has a large isotropy group. In simple cases it is enough to consider 
the length of the orbit of x, namely #{gx : g £ G} to distinguish different isotropy 
types. We now discuss the isotropy groups that occur for the action of on labelled 
shapes. 

The shapes with the highest symmetry are regular pentagons, which are the penta- 
grams = ±(f , f , f , f , f )* and the regular convex pentagons = ±(f , f , f , f , f )*• 
These special shapes are fixed under the whole group , that is their isotropy group 
is -D^, and their orbit length is 1. 

The next group of symmetric labelled shapes are reflection symmetric, with respect 
to some axis through the centre of mass, which is fixed at the origin. After rela- 
belling the vertices by some power of R the relative angles of a symmetric shape are 
{ipi, ip2, ip3, ip2, fpi), which is clearly fixed under MV. Hence the isotropy group of these 
shapes is Z2 generated by MV (or R^MVR~^), and their orbit length is 5. Examples 
of reflection symmetric shapes are shown in Figure [9] at t = 0,r/10,r. 

All other labelled shapes have a trivial isotropy group, so their orbit under 
is length 10, and they do not posses any symmetry. 

The fundamental region of the action of is constructed in the a-space by using 
the reflection symmetric shapes as a boundary. The action of V and M on (01,02,03) is 
simple: y(ai, 02, 03) = (— ai, —02, 03) and M{ai, 02, 03) = (— ai, —02, —0:3) = (ivr — 
ai, ±71 — 02, '^ — Ci3) where the it sign is negative if Oi is negative, and positive otherwise. 



then by using the equivalence relation (11) we get MV{ai,a2,a3) = (01,02,— 03 



(ivr -|- oi, ibvr -|- 02, vr — 03). Thus labelled shapes with 03 = are fixed under MV. The 
operation R in explicit terms is somewhat more complicated, but we do not require that 
formula. 
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Figure 3: The basic symmetry curve 03 = that is used to generate the division of 
the labelled shape space. One quarter of it, denoted by 61, can be transformed into the 
curve 63 by the symmetry operation V 01 M and vice versa; similarly with the curves 
62 and 64. Notice how all the bi curves are ^ rotations of each other. The superscript it 
signs on the pentagons indicate orientation. 



Note that 03 = implies ^2 = V'4 and "01 = V'S; which is an expression of the 
reflection symmetry with respect to vertex 3. Each point on the curve 03 = is fixed 
under MV, while M and V leave the curve invariant as a whole. In fact the action of 
M and V on the curve 03 = are the same, simply (ai, 02) (— ai, —02)- We denote 
the curve 03 = and its images under as symmetry curves. These curves are fixed 
sets of involutive elements in D'^ of the form R^MVR~'^, where k = corresponds to 
the symmetric shapes for which vertex 2:3 is on the symmetry line, since MV fixes "03 • 

The orbit of refiection symmetric shapes under has length 5, and hence there 
are 5 symmetry curves obtained by letting R the generator of -D^/Z2 act on the basic 
symmetry curve 03 = 0. These 5 curves divide the surface of labelled shapes into 10 
simply connected regions, since 10 is the length of the orbit under of a generic (that 
is non-symmetric) labelled shape. 

Explicit formulas for the basic symmetry curve a3 = are obtained from C(ai, 02, 0) = 
0. Figure [3] shows a quarter of the basic symmetry curve parametrized as 



bi{t) = ( t, arccos 



^(3 + 4cost)sec^ (I 
o \ z 



, I where t 



2tt 47r\ , , 



and the subscript a is used to denote the a-coordinate system. 

When t = the relative angles are (-f , f , 0)„ = (f , f , f , f , f )*, this gives 



the positively oriented pentagramj^ When t = the relative angles are (^, ^,0)^ 
^(0:1,(12,0:3)0, ~ (^1, i/'2, V'3i V'4i ^5)' shows equivalent relative angles in different coordinate systems. 
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10 1 

2 2 



0(1 

Figure 4: Projection of the shape space onto the plane (ai, 02) with contours of constant 
moment of inertia. Taking the rotation symmetry R'^ for each of the blue curves bi 
from Figure [3] produces the red, green, yehow, black curves denoted by rj, gi, in, bi, 
respectively, for i = 1, 2, 3, 4. These are the fixed sets of R^MVR'^ for A; = 0, . . . , 4 with 
vertex z^-k fixed when defining the vertex reflection. Of the 10 topologically equivalent 
regions we choose the piece bounded by the curves r4,g3, k2 and yi as our fundamental 
region and denote it by cj). 



(x'X'X'X'T")*' ^^^^ corresponds to the positively oriented regular convex pentagon. 
Applying M oi V to bi gives 63. In the plane 03 = this amounts is a rotation by vr. If 
instead a rotation by ^ is performed the curves 62 and 64 are obtained. This is a result 
of the fact that the equation 0(01,02,03) = in (13) is even in Oi. 

The images of the basic symmetry curve of shapes for which vertex Z3 is on the line 
of reflection symmetry in the plane are shown in Figure |4j The basic symmetry curves 
are the fixed set of the involution MV. The /c*^ image of this curve under R is the fixed 
set of R^MVR~^ , with vertex z^-k on the line of reflection symmetry in the plane. The 
colour scheme for vertices and symmetry curves was chosen, such that shapes on the 
symmetry curve with colour c, have the vertex with that colour on the line of reflection 
symmetry in the plane. 

We chose a fundamental region as shown in Figure |4| which is bounded by R{b4), 
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Figure 5: Schematic picture of how to glue the labehed shape space from copies of the 
fundamental region (p. The labelled outer edges with their specified orientation are to be 
glued together to form the shape space of genus 4. The symmetry operations required 
to map (p to any other region is identified in this diagram. 



i?^(62)j and R^^b^). The fundamental region cp has the useful property that it can 
be parameterised by ai and 02- Starting with (j), the nine remaining regions are obtained 
by applying the symmetries R and MV to (p. When considering an arbitrary point in 
shape space, it may always be moved to the fundamental region (p by applying some 
symmetry operation from . The topology of how the 10 pieces are glued together to 
give all of the shape space is shown in Figure [5} This figure is somewhat reminiscent 
of a tiling of the hyperbolic plain obtained from the triangle group of type (2, 5, 5), but 
since the tiles are quadrilaterals in our case the analogy is rather incomplete. Now all 
the analysis required can be done in (p. The remaining parts of the surface are covered 
by the orbit of (p. 

In Figure [4] (and later figures) two copies of each orientation of the regular pentagon 
and pentagram are shown. Note that these are equivalent under (11), so they represent 



the same point of the labelled shape space. The reason for displaying both is to emphasise 
the peculiarity of our fundamental region, which occurs because the corner points of the 



4 MOMENT OF INERTIA AND ANGULAR MOMENTUM 



14 



fundamental region (j) are actually outside of the fundamental domain of the a's since 
as G [0,7r). 

Finally, we discuss the reduction by the remaining symmetry V, which in general 
reverses the orientation of a given oriented shape. Our fundamental region (j) contains 
both of these, since they are different as oriented shapes. The quotient by V identifies 
these oriented shapes with different orientations as the same, and the result is the set 
of un-oriented shapes, or just shapes in standard terminology. A natural way to define 
a fundamental region that contains each (un-oriented) shape once (up to congruence 
E{2)) is to cut (j) into halves along ai = a2- Along this line the total angle sum Yl'^Pi 
jumps from vr to — vr. We denote (/>"*" as the half of (p that contains each shape with a 
positive angle sum. Thus 0"*" contains each (un-oriented) shape exactly once. Half of 
the additional side ai = a2 of 0"*" is open while the other half, including the midpoint, 
is closed. Every point in the interior of 0"*" is generic and has a Diq orbit of length 
20. The origin has orbit length 10, while the sides have length 10 and the corners have 
length 2. Hence there is one new symmetric shape, which is invariant under V (but not 
under M), corresponding to the origin ai = a2 = of 0. This peculiar shape is shown 
in Figure 6(c) and is invariant under orientation reversal. 



4 Moment of Inertia and Angular Momentum 

The moment of inertia of the equilateral pentagon with point masses mj = 1 at the 
vertices with respect to its centre of mass is given by 

1=1 

which is obviously invariant under the symmetry group -Diq. In the a-coordinate system 



equation (16) becomes 



8 12 6 

1 = 4 + 2 cos ai cos 02 H — cos ai cos 03 H cos 02 cos 03 H — sin ai sin 02 ■ (17) 

5 5 5 

The contours of constant / are shown in Figure |4} Notice how ([Tt]) is even in 03; this 
implies 03 can be eliminated using (14) to give an expression that is a rational function 
of trigonometric functions. The denominator of this function of two variables vanishes 
at the 8 points given by {±tt/6, ±57r/6)a and (ib57r/6, ±7r/6)a where taking all possible 



combinations of the it sign; this has already been discussed in relation to (14). These 
points are all outside the fundamental region so they do not cause a problem. The 
moment of inertia is invariant under the full symmetry group Diq. To find all critical 
points of the moment of inertia we employ symmetry reduction, and hence considered 
only the critical points of the moment of inertia within the fundamental region . The 
action of Diq then generates all critical points on the entire surface. On <j), as shown 
in Figure |4j there are five critical points of the moment of inertia. Four of them are on 
the boundary, and are in fact on the four corners and one is in the interior. The four 
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(a) 2 minima; 
J= i(5- V5); 
isotropy: {R, MV) 



(b) 2 maxima; 
/= 1(5 + 75); 
isotropy: {R,MV) 



(c) 5 saddles; 

^ 2 ' 

isotropy: (V) 



(d) 5 saddles; 
I = ^■ 

isotropy: {V} 



Figure 6: Applying the symmetry operations R, M and V to the critical shapes shows 
that there are in total 14 critical labelled shapes. The types are: 2 minima [6(a) 



maxima [6(b)], 10 saddles (5 of the form |6(c)[ and their mirror images under M 6(d)|) 



critical points on the boundary of the fundamental region also happen to be located on 
the boundary 03 = vr where (01,02) are not local coordinates. Thus the criticality is 



established using (17) and Lagrange multipliers incorporating the constraint (13) 



2(5 — \/5) and are located at ib(- 



The pentagrams in Figure 6(a) have global minimal moments of inertia with value 



and are located at ib(|^. 
shown in Figure 6(c) and 



■ 5 ' 



vr 



_ _ _ TT 7r^ 

5 ' 5 ' 5 ' 5 ' 5 ' 



The regular convex 



pentagons in Figure 6(b) have global maximal moments of inertia equal to 2(5 + ^/b) 

^,Tr)a = ±{^,^,^,^,^Y. The saddle critical shapes 



6(d) have moments of inertia equal to |; and the saddle found 



in (f) occurs at (0, 0, vr 



vr 



k), K, vr, — K, — 2(vr — k)) where k = arccos(|). 



Applying the symmetry operations to this saddle produces a total of 10 saddles on the 
entire surface. Recalling the Dio orbit lengths of the symmetric shapes there are 2 
minima, 10 saddle points, and 2 maxima, so that the Euler characteristic of the shape 
space isx = 2 — 10 + 2 = —6 proving that the genus is 4, as claimed earlier. 

For a system of particles with coordinates (xj, yi) in the plane the angular momentum 
with respect to the origin (which is the centre of mass in our case) is a scalar given by 



L = '^mi{xiyi - Xiyi) 



(18) 



where the dot denotes the time-derivative, so that Xi is the x-component of the velocity 
vector of vertex Zi. In the case of the equilateral pentagon where it is placed on the 



complex plane with each vertex having unit mass, equation (18) simplifies to 
5 

^9(ziii) wheie Zi = Xi + iyi . (19) 



L 



i=l 



Computing the angular momentum by transforming the vertices found from equation 
([3]) into the a-coordinates and using equation (19) gives 



L = ie + Fidi + F2d2 + F3Q3, 



(20) 
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where / is the moment of inertia as found in equation (17) and 



~ 4 6 

Fi = 2 + cos ai cos a2 H — cos ai cos 03 H — cos 02 cos Q3 

5 5 

3 6 , , 

H — smaisma2H — sma2sma3, (21a) 
5 5 

~ 4 6 

F2 = 2 + cos «! COS 02 H — COS ai cos 03 H — cos 02 cos Q3 

5 5 

3 4 

+ - sinoi sina2 + ^ sinai sina3, (21b) 

~ 12 8 4 6 

F3 = ^ — cos «! COS a2 H — cos ai cos 03 H — cos 02 cos 03 

5 5 5 5 

8 , , 

H — sin Qi sin 02- (21c) 
5 

The same expression for the angular momentum L in terms of orientation and shape 
coordinates (6,01,02,0(3) can be derived by starting with the Lagrangian given by the 
kinetic energy £ = ^^mi\zi\'^ . After introducing orientation and shape coordinates 
the Lagrangian becomes independent of 9, and the conjugate momentum dC/dO is the 
conserved angular momentum. In this derivation we would either treat 03 as a known 
function given by (14) or use Lagrange multipliers. In the free motion of the system 
the shape and orientation would be determined by the corresponding Euler-Lagrange 
equations. In the following we instead consider the shape as given by explicitly time- 



dependent functions ai{t), and the only equation of motion we use is (20) to find 0{t). 

5 Geometric Phase 

In the coordinates [0, ai, 02, 03) we can interpret the formula for the angular momentum 
as being decomposed into a contribution that originates from changing the shape, and 
a single term 19 that originates from the rotation of the shape. For a rigid body this 
would be the only term present. The well known but nevertheless surprising result is 
that even when L = the orientation 9 may change. This is most clearly seen when we 
solve for 6, which gives 

L I F F F \ L 

^ = J - -jdi + -jd2 + -jds = J + + ^2^2 + Fsds- (22) 

This gives the decomposition of the change of 9 into a dynamic phase proportional to 
the constant angular momentum L, and a geometric phase proportional to the shape 
change, which is proportional to the time derivatives of the angles aj. Henceforth we 
will set the angular momentum to zero, L = 0, so that we can study the question of how 
to maximise orientation change in the absence of angular momentum. 

The geometric definition of the rotation angle 6 is not unique, and making a particular 
choice is called a gauge. Given a particular shape change ai{t) the resulting overall 
rotation 9{t) at a particular time t depends on the gauge. To get a gauge invariant 
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quantity we consider A6 = 0{ti) — 6{to) for closed loops 7 in shape space, that is for 
o:i{ti) = ai(to). For a closed loop in shape space it makes sense to subtract ^'s to get A^, 
since they are computed for the same shape, see for example [7] for more information 
about gauge invar iance. 



Integrating equation (22) with L = over a loop 7 in shape space gives the overall 
change in 9, 



A9= (j) Fidai + F2da2 + Fada^. (23) 
7 

It is not obvious how to choose the loop 7 such that this integral becomes large. In order 



to find the optimal loop we use Stokes' theorem to change the line integral (23) into a 
surface integral 



= yy V X F • dS (24) 

Si-y) 

where F = (^1,^2,^3)*, Vx is the curl with respect to the angles Oi, and S{'j) is the 
surface enclosed by 7 on shape space C = 0. The formula for A9 is an integral over a 
two- form (after eliminating, say, 03). We want to convert this integral over a two- form 
into an integral over a function on shape space. The function thus defined will be called 
the magnetic field B. In order to do this we need a metric g on our shape space surface, 
so that the magnetic field becomes invariantly defined after dividing the two-form by 
the area form ^/det g. 

In general a "magnetic field" arises from reduction by a continuous symmetry. Like 
the equation for the angular momentum L, the magnetic field could be derived in the 
Lagrangian or Hamiltonian formalism, see for example |9l [7]. There, one would use 
e.g. (ai,a2) as local coordinates and obtain a magnetic field (V x F) • VC/C3 where 



C3 = dC/das is the derivative of the constraint (13). The natural metric g in our case 



is obtained from the kinetic energy of the Lagrangian of the Pentagon after symmetry 
reduction, so that = (V x F) • V C / {C3^/detg) is the magnetic field. However, the 
expression for detg is fairly complicated, and since we are not interested in the free 
motion of the pentagon (which is given by the geodesic fiow of this metric) we prefer to 
use a different approach. Suffice it to say that the magnetic field B is invariant under 
our discrete symmetry group Diq and has zero average over the whole shape space. The 
geometric interpretation of this average as a Chern class confirms that the bundle of 
orientation 9 over shape space is a trivial bundle. 

Instead we attempt to use the metric induced by the embedding of shape space 
in the three-dimensional torus T'^ with coordinates Oj given by the equation C = 0, 
see (|13l). The area form of this surface using qi and a2 as local coordinates is dS = 



IVCI/C3 daida2, and with dS = ndS where n = VC/|VC| the integrand in (24) becomes 
(V X F) • VC/C3 daida2 = (V x F) • VC/\VC\dS, so that the magnetic field would be 
the factor multiplying dS. This function differs from i? by a scalar non-zero factor. The 
trouble with the function (V x F) • VC/|VC| is that it is not invariant under the discrete 
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symmetry group. This is not too surprising, since we arbitrarily invented a metric on 
shape space by embedding it into T'^ using C = 0. In the following all we want to do 
with the magnetic field is to study its zero level curve B = in the fundamental domain 
(f). Thus multiplying i? by a non-zero factor does not make a difference. In order to 
make the result invariant under the discrete symmetry group we multiply the function 
obtained from the metric of the surface C = by |VC|, and thus define our magnetic 
field to be given by the simple expression 

B = {V xF)-VC. (25) 

Explicitly we find 
5 

-I^B = cos(3a3) [ cos qi + cos 02] + cos(2a3) [3 + 2 cos(ai — 02) + cos(ai + a2)] + 
8 

+ cos as [ cos ai + cos 02 — cos(3q;2) — cos(2ai — 02) — 2 cos(ai + 202)] + 
+ sinai sin 02 — cos ai (2 cos 02 + cos(3a2)) — cos(2ai) — 2cos(2a2) ■ 

The "magnetic field" B is invariant under the the action of the discrete symmetry group 
DiQ. We don't have a good explanation why this expression possesses the correct discrete 
symmetry, but this together with its simplicity if the main reason to use it instead of 
B. Using the proper magnetic field B instead of B does not change our results, only a 
different parametrisation of the same loop B = would be used. 



6 Optimal Shape Change 

We are seeking the optimal contractible loop 7 on shape space in the sense that it 



maximises overall orientation change of the pentagon given by ( 23 ) after one revolution. 
Imagine we start with a small loop 7 in a region where B > 0, say near a positive 
maximum of B. Enlarging the loop will increase \ A6\ as long as -B > 0. This process can 
be repeated and we can keep growing the loop, yielding a larger and larger \A9\ . However, 
enlarging the loop 7 across the B = contour would give an opposing contribution 
towards the integral, thus lowering \A9\. Hence the largest \A6\ is achieved when the 
loop 7 coincides with the contour B = 0, assuming there is a contractible zero-contour 
of B enclosing the initially chosen small loop. In general, the connected components of 
the B = contour on a genus 4 surface may be non-contractible. Yet in the present case 
all zero-contours of B are contractible. 

The sign of is dependent on the orientation of 7. A positively oriented loop in 
the B > region will yield a negative AO. Though this sounds counter intuitive, the 



reason is that the matrix A from the affine transformation given in (10) has a negative 
determinant. That is, a positively oriented loop in a-space corresponds to a negatively 
oriented loop in -space. The construction works similarly starting in an initial region 
with B < 0, where a positively oriented loop will yield a positive AO. 

The connected components of the zero-contours of B give analytic curves that yield 
locally optimal loops that make AO extremal, in the sense that any small variation of the 
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Figure 7: The light blue region indicates B > and pink B < 0. There are 4 curves with 
B = in the fundamental region, enclosing the positively/negatively oriented regular 
pentagon/pentagram, respectively. 

An interactive Mathematica file link to manipulate, cdf allows you to select a point in the 
fundamental region (j) with the mouse cursor and displays the corresponding pentagon. 



loop decreases the value of |A^|. Globally it is possible to connect two disjoint regions 
with B > with a curve through the area with B < 0. Traversing this curve back and 
forth provides no net contribution as it does not enclose any area. Alternatively one 
could also traverse the same loop twice and hence double the amount of rotation. 

It is enough to consider the B = contour in the fundamental region (p because of 
the discrete symmetry, as shown in Figure [7| The symmetry V maps curves enclosing 
the positively oriented regular pentagon/pentagram into curves enclosing the negatively 
oriented regular pentagon/pentagram, respectively. As V preserves A9, only the two 
loops shown in Figure [7| need to be considered. The orientation is as indicated and 
specifically chosen such that the overall rotation A9 is positive. 

When defining the optimal loop 7 through B = we may multiply B by any smooth 
function on the shape space that is invariant under the symmetry group Dio without 
changing the result. The overall denominator of B is the moment of inertia squared. 



which can be removed when computing 7. The integrand in (24) involves rational func- 



tions of trigonometric functions of Oj that is even in a^. Thus we can replace 03 by 
03(01, 02) and obtain a rational function of trigonometric functions of ai,a2- We have 
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Figure 8: The diagram shows how the relative angles ipi change along the B = contour 
around the pentagram, the colouring of V'i is from our standard colour code, that is, 
i = 1, . . . , 5 is green, red, blue, black, yellow. The curve segment in (p as seen in Figure 
[7] was used to construct the entire loop on shape space using the indicated symmetry 
operations. 



not been able to find a simple closed form for this integral, so we resort to numerical 
methods at this point. 

Once the curve is known in the fundamental region the symmetries are used to find 
the curve in all of shape space. The result of this is shown in Figure [Sj In order to 
parameterise the zero-contour of B, we numerically solve Hamilton's equations with a 
Hamiltonian given hy H = BI^ with initial conditions on the boundary of the funda- 
mental region (p. Note that this Hamiltonian is merely used for the purposes of obtaining 
a parametrisation of the B = contour; it is not the Hamiltonian of the free motion of 
the pentagon. The speed with which the contour is traversed is irrelevant for the final 
geometric phase A9, so whether we consider H = B or H = BI^ makes no difference. 
We can even multiply by -B by a function that is not invariant under the Diq action, but 
then the overall solution pieced together from the action of the symmetry group on the 
fundamental piece may not be smooth but only once differentiable. The ODE in local 
coordinates (01,02) is 

dH . dH 

ai = -b- — , a2 = b- — . 
oao om 



Here H = H{ai, 02,03(01,02)) and S = C3(ai, 02, "3(0!) 02)) is the symplectic mul- 
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tiplier from the area-element. Maybe the easiest way to derive this non-standard sym- 
plectic structure is to start with the global ODE on the three-torus a = VC x Vi? 
(which also gives an alternative way to numerically compute the whole loop B = 0) and 
consider VCx as a Poisson structure. Reduction to the symplectic leaf {C = 0} using 
local coordinates (01,02) then gives the symplectic structure Sdoi A da2- 

The solution curves Oi{t), and hence ipi{t) as well, along the B = contour are 
periodic with period r. The solutions for the different angles are related by a phase 
shift, thus it is enough to study a single curve for the whole period r, say ip5{t), which 
is an even function. Specifically the positively oriented loop in a-space around the 
positive pentagram has the relation ipi{t) = ipi-^-i{t + |) while the negatively oriented 
loop has relation ipi{t) = ^i+i(t — |)- The phase shift relation for the positively and 
negatively oriented loops in a-space around the positive regular convex pentagon are 
iJiit) = ipi+i{t + ^) and i{ji{t) = ipi+iit - x)' respectively. 

A natural way to encode the final answer is to decompose ip5{t) into a Fourier cosine 
series, 

°° 27r 
V'5(*) = / a„ COS (nwt) where a; = — = 7.3634..., (26) 

n=0 

and the coefficients a„ are given in Table [T] 



n 


an 


n 


an 





0.6283 






1 


-0.9646 


8 


-0.004078 


2 


-0.3974 


9 


0.0003710 


3 


0.1595 


10 





4 


-0.0779 


11 


-0.0002426 


5 





12 


0.0001005 


6 


0.008748 


13 


0.00003240 


7 


-0.002677 


14 


0.00001171 



Table 1: Numerical values of the Fourier coefficients On of ip5{t). Along this loop where 
-6 = truncating the Fourier series at n = 14 produces an absolute error \B{ipi{t))\ on 
the order of 10^^. 

The fact that all coefficients a^n = is equivalent to J2^=o ^« + ) — const = 
5do. Using the property ipi{t) = ^l'^+l{t + |) translates this into Yl'^Pii'^) — ^^o- By 
construction we have that '^ipi = it mod 2tt, hence do = 7r/5. If we consider the sub- 
Fourier series defined by fi{t) = ai+^n cos((Z+5n)a;t) so that /o(t) = |, we find that all 
angles can be explicitly written as linear combinations of four functions fi, / = 1,2,3,4 
(and the constant function /o) as 

4 

V^5+i(t) = 5ii;//(i)e'^^^''/^ 
1=0 
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The function fi{t) can be written as e^'^^^ fi{t) where fi{t) = X] Q+fcnC*^'^"* has period |. 
In real variables this becomes fi{t) = cos{lujt){aifi+ai^i cos(5a;t)+. . . ) where ai^n = a«+5n- 
Any periodic shape change that is obtained by unfolding a curve from the fundamental 
region has these properties. 

As ^ is a composition of periodic functions, it is also a periodic function with the 
same period r, with Fourier series 

oo 

9{t) = ao + {a„ cos [nuji) + bn sin (nwt)}, (27) 

n=l 

where the numerical values of a„ and bn are given in Table [2} Notice that even for the 
optimal loop the sign of 9 is not constant along the loop. 



n 


an 


bn 





0.9239 




1 


-1.8671 


-5.7463 


2 


-2.4891 


1.8085 


3 


-1.4990 


-1.0891 


4 


0.6032 


-1.8566 


5 


-0.3977 





6 


0.1016 


0.3127 


7 


-0.0587 


0.04265 


8 


0.1022 


0.07423 


9 


-0.006463 


0.01989 


10 


0.005644 






Table 2: The Numerical Values of the a„ and b„ terms of 9 



Integrating equation (27) term by term gives 



oo ^ 

9it) = C + aot + ^l 

n=l ^ 



sin (nut) — cos (nujt) \. (28) 



nu) nuj 



Hence = 9{t) — 0(0) = a^T. The numerical value for obtained along the locally 
optimal loop around the regular pentagram is 

A6' = aor 0.78837 w 45.17°. (29) 

Performing the same analysis for the locally optimal loop around the regular convex 
pentagon yields 

A6' pa 0.49147 28.16°. (30) 

We conclude this section by illustrating how the shape of the pentagon changes as 
the loop of optimal shape change is traversed, see Figure [9| Although it may appear 
that the yellow vertex travels along the blue-red edge this is only approximately true. 
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(g)t=ij^ (h)t=^ (i)i=^ ii)t=m (k)i=T§5 mt-r 

Figure 9: The optimal shape change of the equilateral pentagon that maximises overall 
rotation at angular momentum zero after one period. The diagram shows the shape 
change within the fundamental region (j) only (except for the last image which serves for 
comparison between the initial and final orientations of the pentagon). By symmetry 
reconstruction, the remaining shape changes outside the region continues on with this 
motion in the reverse order; then the initial motion shown is repeated, each time with 
different colouring of the vertices until it returns to its original state. 
Movie files links and/or the Mathematica link to animations, cdf file shows the complete 
animations along the different optimal loops of the equilateral pentagon. 
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7 Convex Pentagons 

The results obtained so far apply to the family of all equilateral pentagons including 
degenerate and non-simple pentagons. We now restrict to the sub-manifold of equilateral 
convex pentagons. The intuitive approach to finding the sub-manifold containing all 
equilateral convex pentagons from the original manifold is to consider the boundary 
cases for which the pentagon is not strictly convex. Fixing one of the angles of the 
equilateral pentagon to vr gives the curve on shape space which contains all possible 
equilateral pentagons with the chosen angle ipk = for some A; G Z5. The importance of 
this curve is that parts of it separate the strictly convex pentagons from the nonconvex 
pentagons locally. Choosing '03 = tt we can see from ([To|) that ai = 02, so one way to 



parametrise this curve is to set 03 = —t, substitute in ( 13 ) and solving for ai = 02 gives 



p{t) = {arccos xit) , slvccos x{t) , —t} where x(0 = — cos t + y cos^ t — -. (31) 
To orient the curve in the clockwise direction as t increase, we included a minus sign in 



the third component of the parametrisation in (31), the first two components are even 
functions so the minus sign may be omitted. 



To find the endpoints of the interval where (31 ) gives the boundary separating convex 
and non-convex pentagons, we need another angle ^|Jj, j 7^ 3 to equal vr. This occurs when 

t = ± arccos (^|^) = ^Ps, thus t G [—Ps,Ps] is the interval of interest. The pentagonal 
shapes at these endpoints are the isosceles triangles with side lengths {1,2,2}. 

Note that when t G [—Ps,Ps] the boundary curve traverses through two different 
fundamental regions, surrounded by the symmetry curves {^1, ^2, fcs, 1/4} for t G [—ps,0] 
and {bi, k2,r3, g^} for t G [0,^^]. These two fundamental regions are related by the 
symmetry operations MV. To complete the white loop as shown in Figure [10(b)] four 
additional symmetry operations i?*, i G {1,2,3,4} need to be applied to p{t). In order 
to keep the curve in a single fundamental region we restrict t to [— p^, 0], so that R^{p(t)) 
for t G [— ps,0] gives the black curve inside (j) with B > as shown in Figure [IT} The 
shape at i = is the trapezium with side lengths {1,1,1,2}. 

The orbit of p{t) under gives the piecewise smooth loop around the regular convex 
pentagon shown in white in Figure |10(b)[ The particular order of symmetry operations 
can be read off from Figure [5j The corner points where the loop is non-smooth are the 
isosceles triangles that have two angles ipj equal to vr, while the midpoint of the white 
edges indicated by the coloured points correspond to the trapeziums with side lengths 
{1,1,1,2}. The corresponding loop enclosing convex pentagons of the opposite orientation 



is obtained by applying either M or V. Also shown in Figure 10 are the B = contours in 
black. Intersections of the loops with the boundaries of the 10 fundamental regions that 
meet at the central positively oriented pentagram (left) or positively oriented regular 
convex pentagon (right) are marked by coloured dots, also compare Figure [5] 

The white loop which demarcates the boundary of the region of convex pentagons lies 
completely inside the region with B > (given by the black loop). Thus the maximal 
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(a) 



(b) 



Figure 10: Two views of shape space C = showing the complete optimal loops. The 
gray point in the interior of the loops denotes the pentagram (left) and regular convex 
pentagon (right). Each loop intersects the boundaries of copies of the fundamental region 
at 10 different points denoted by the coloured points, the colour is chosen to match the 
corresponding symmetry curve, compare Figure [5j The dashed white lines denotes the 
intersection with a,- = ztyr. 



rotation of the equilateral convex pentagons is obtained along this loop. The argument 
is similar to that used to show optimality for the B = loop. Deforming the loop into 
the non-convex region is not allowed, while deforming the loop in the other direction 
would decrease the enclosed area, and hence the overall amount of rotation as > 0. 

The derivatives ipi along the boundary of the convex region have finite jumps at kps 
when k is odd, these jumps are shown in Figure |12(b)| and corresponds to the corners 
of the white loop from Figure 10(b) , The corresponding shapes at these corners are 



the isosceles triangles with side lengths {1,2,2}. For even k in kps, the shapes are the 
trapeziums with side lengths {1, 1, 1, 2} and the transit ion fro m one fundamental region 



to another at this point is smooth, as seen in Figure 12(b) By integrating 9 in (22) 
we obtain a continuous function 6(t) (though not smooth at kps for odd k), and the 
overall rotation generated by this loop is found to be A9 ~ 0.33117 ~ 19.01°. This is 
the maximal overall rotation achievable after one period if we restrict the allowed shapes 
to the space of all equilateral convex pentagons. 
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The optimal way for the equilateral pentagon at zero total angular momentum to achieve 
maximal overall rotation after one period of a periodic contractible loop is to follow the 
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4 



Figure 11: This is an enlargement 
of Figure [7| In addit ion the oriented 
black curve inside the blue region 
indicates the boundary of convexity 
given by R^{p{t)). The symmetry 
operator is applied to p{t) to 
map the curve to the fundamental 
region 4>. 



ai 




Figure 12: The relative angles ipi (left) and (right) along the loop of the boundary 
of the convex pentagons in the clockwise direction starting at the red point located at 
p{—Ps) in Figure [lO(b)| where p is from ( |31[ ). The colouring of the curves i = 1, . . . , 5 is 
green, red, blue, black, yellow (standard colour code). The dashed vertical lines denotes 
the boundary of each fundamental piece, and the finite jumps for ipi that occur at every 



second boundary are the result of the corner points from Figure 10(b) 
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B = contour around the regular pentagram. This results in an overall rotation of 
AO M 0.78837 radians or 45.17°. 

If we restrict to the subset of simple equilateral pentagons, then the B = contour 
around the regular convex pentagon provides the maximal overall rotation with A9 ~ 
0.49147 radians or ~ 28.16°, this is the black loop found in Figure 10(b) Finally, if we 
further restrict to the subset of equilateral convex pentagons, then the maximal overall 
rotation reduces to A9 ~ 0.33117 radians or ~ 19.01°. The fact that this loop is smaller 
and entirely contained within the larger B = loop around the regular convex pentagon 
(see the white loop in Figure |10(b)[ ) explains why the overall rotation is smaller for this 
loop. 

An intuitive explanation of why the loop around the pentagram gives a largest value 
is that the moment of inertia for the pentagram is the global minimum. Hence the 
magnetic field B tends to be bigger there, as compared to the magnetic field near the 
regular convex pentagon, which has maximal moment of inertia. 
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Appendix 

Movies of .mp4 format for optimal shape changes are available at 



http : / /www . maths . usyd . edu . au/u/ williaiiit/peiitagon/BeqO_peiitagrain . mp4 



for the optimal loop around the pentagram, 



http : / /www . maths . usyd . edu . au/u/ williamt/peiitagon/BeqO_convex_peiitagoii . mp4 



for the optimal loop around the pentagon, and 



http : / /www . maths . usyd . edu . au/u/ williamt /pentagon/subset _convex_pentagoii . mp4 



for the optimal loop for convex pentagons. 
Interactive Mathematica CDF files are available at 



http: //www. maths .usyd. edu. au/u/williamt/pentagon/manipulate . cdf 



for finding the shape corresponding to a point in the fundamental region, and 



http: //www. maths .usyd.edu. au/u/williamt/pentagon/animations . cdf 



for animations of the optimal shape changes. 
Download the Wolfram CDF player at 



http : //www . wolfram . com/cdf -player/ 
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